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In this work critical speeds of rotating shafts 
are obtained. Isotropic and anisotropic (dissimilar 
stiffness in two principal directions) shafts are onsidered. 
Various factors affecting the critical speeds are taken into 
account and their effects are studied. 

The relevant governing differential equations are 
derived and FEM has been used for their solution. G-alerkitf's 
technique is used to minimise the residues obtained by finite 
element approximation. The element matrices are integrated 
numerically using Gauss quadrature, assembled and solved for 
frequencies and critical speeds for various boundary condition 

The factors affecting the critic'- 1 speeds considered 
for the present study are shear deformation, rotatory inertia 
and gyroscopic moments. Shaft with dissimilar support 
conditions in two principal planes has also been studied. 

:.n attempt to estimate the errors due to the assumptions of 
massless shaft and rigid disk has been dene. The critical 
speeds for multiple span shafts are also obtained. 

The effect of internal damping on the stability 
has been studied. Both viscous form and hysteretic form of 
internal damping are considered. The effect of bearing chara- 
cteristics (stiffness and damping) on the instability threshol 
has been studied. 


CHAPTER- I 


INTRODUCTION 


1 .1 THE PROBLEM IN BRIEF : 

It was recognised early in 18th century that, 
a rotating shaft initiated vibrations of large amplitudes 
at certain speeds and often caused catastrophic failure of 
supporting structure and other parts of the system. The 
operating speeds at which this phenomenon occured, came to 
be knowm as critical speed. 

A’ shaft always has certain amount ol mass unbalance, 
which during rotation, deflects the shaft axis from its 
eeuillibrium position. This leads to trsnsverse vibration 
of shaft at its natural frequencies. In general the natural 
frequencies of a rotating shaft are dependent on speed of 
rotation. The vibrations of shaft along any two mutually 
perpendicular directions in the plane of cross section, 
occur at same frequency, but with a phase difference of + 

Due to this ? the shaft axis begins to precess about the line 
joining the supports, with a speed equal to the angular 
frequency of natural vibrations. The precession of shaft 
can be in the same direction as that of spin of shaft or in 
the opposite direction. V/hen precession is in same direction 



£.2 the t of spin, it is tErmed as 


forward precession and 
viiea in opposite direction, it is called as backward precession 

The inertia force due to unbalance has the frequency 
of spin speed of shaft. Hence, when the speed of precession 
of shaft becomes equal to the spin speed, resonant state is 
reached, i.e. amplitude grows linearly with time. When this 
state arises, the corresponding spin speed of shaft is 
called as critical speed, Depending upon the direction of 
precession, it is termed as forward critical speed or 
backward critical speed, 

1 .2 LITERATURE; REVIEW : 

The first published work related to this problem 
was in 1669, by Rankine [1 ], The operating speeds of rotating 

machinery were below the first critical speed in those days . 
Hence, most of the early work in this field is limited to the 
determination of first critical speed. However, with the 
advent of new technology in the field of rotating machinery, 
the present d^y operating speeds are f:r in excess of first 
critic: 1 speed. Consequently, the more recent literature 
treats a greater range of problems related to rotating shafts,, 
such as effect of various system parameters on critical 
speeds, determine tion of stable operating speed ranges, 
causes of instabilities and methods to improve the stabilities . 
Thus the study of vibratory motion becomes necessary, for 
stable and safe operation of rotating machinery. 
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a very good survey of work done in this field till 
1967, can be found in the monograph "Dynamics of rotating 
shaft” by Loewy and Piarulli £2 J. This monograph tries to 
illustrate the basics of rotor dynamics . It neglects the 
mass of the shaft, i .e ., it uses lumped mass approach. It providi 
general but qualitative information on the effects of 
v rious parameters on critical speeds 'and stability. It. 
also treats coupled, lateral and torsional motion and 
balancing . 

1 more involved and analytical study of the problems 
of rotor dynamics is found in the work of Tondl [3 ] . The 

book gives c good insight of several phenomena of rotor 
dynamics. The topics emphasised in his work are, effect of 
internal damping, effect of unequal shaft stiffness and effect 
of. oil film in journal bearings. It deals comprehensively with 
problems concerned with the stability of rotor motion, with 
the creation of self excited vibration, and _with non-linear 
resonances. But almost the whole book is based on the 
assumption of massless shaft and rigid disk. All the same, the 
book gives an excellent review of the problems, supplemented 
with results of experimental investigations. Combined effect 
of torsional and flexural rigidities, combined effect of 
several factors, causing instability have been dealt in 
detail . 
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Another very prominent contribution in this field 
is of bimentberg £ 4 ] , This worl: denis with flexural 

transverse vibration, Major pert of the book uses continuous 
approach, closed form solutions, for hinged -hinged support 
conditions, for several cases are obtained. This book showed # 

J 

the limitation of critical speeds in lorvcrd precession of a 
shaft with uniformly distributed mess, was a result of 
neglecting the effect of shear deformation. Effect of 
external and internal damping, effect of dissimilarity of 
stiffness, effect of dissimilar fixing conditions in two planes 
multi -bearing systems, are few of the cases studied in his 
work. Further, both uniform stiffness and dissimilar stiffness 
shafts are analysed with and without consideration to shear 
deformation effect. 


Another recent text in this field is of Rao [5 ] 
It uses lumped mass approach . 


An excellent collection of papers on several 
problems of rotor dynamics, presented at an international 
symposium held at Denmark in 1974, can be found in the book 
edited by Niordson [6 ^ ' 

The application of FEK to the problems related to 
rotor dynamics was done by Ruh3. and Booker [7 ] , in 1972. 

They analysed stability and unbalance response of a turborctor 
system. But the model used, included only translational 
kinetic energy and elastic bending energy i.e. used Euler 
beam theory. Nelson [g ] improved the model by incorporating 
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the effect of rototory inertia end gyroscopic action. 
Formulation is done using energy method. The effects ol bearing 
stiffness and damping are also incorporated . Zorzi and 
nelson [9 ] , incorporated the effect of internal damping 
in the finite element model of [ 6 ] • --ccording to [4 ] 

the possibility of instabilities in the post -critical speed 
range due to internal damping, was first shown by Kimball 
and Niwkirk [ 10 ] * in e chanism of internal damping is 
better represented by considering the hysteretic form, rather t 
viscous form [ 4 ] . 4 combined model incorporating both the 

forms of damping is formulated by [ 9 ] • A study of 

effect of bearing stiffness and damping characteristics on 
the instability threshold is done by them. Later in 1965, 
Greanhill , Bickford end Nelson [ 11 ] , extended the work 

of [ 9 j t for a shaft with variable cross section, by using 
a conical beam finite element. 

There are a good number of authors who use transfer 
matrix method in this field. Lund [12 ]» studied the 

effect of fluid-film journal bearings and hysteretic internal 
damping, on the critical speeds and instability threshold. 

He also incorporated the effect of bearing stiffness and 
bearing damping. But he treats the shaft ss a uniform beam, 
noglecting the effect of shear deformation , rotatory inertia 
and gyroscopic moments. Bansal and Kirk [13 ], improved this 
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\<ork by incorporating the disk gyroscopic moments and 
shear deformation effect. Ramakrishnan and Prabhu [id ] 
have studied the unbalance response, instability and critical 
speeds of multispan rotor systems , 

i-J.1 these autnors, including Rao [ 5 [[use lumped mass ap 
roach in the transfer matrix method. Lumped mass approach is 
* alvays less accurate than the distributed mass approach. 
Moreover this method requires some skill and experience as 
there are chances of missing some modes . Satisfaction of 
boundary conditions of intermediate supports is not very 
straight forward . 

V.’ith continuous approach in transfer matrix method 
for beam problems, one faces some numerical problems. 

Moreover it is quite difficult (if not impossible) to obtain 
transfer matrices for rotors with rotatory, shear and 
gyroscopic effects. This difficulty is easily overcome in 
transfer FEK (TFEM) at the cost of little accuracy, Gupta [ 15 
and Gael [16 ] .It eliminates the need of the solution 
over an element as it uses' the finite element shape functions . 
Further, TFEM also has the advantage of lesser computer 
memory requirement as, the order of matrices to be handled 
remain same as in transfer matrix method . Thus the TFEM blends 
the advantages of FEK and transfer matrix method . But 
obtaining the complex frequencies for damped cases does not 
seem possible. 
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ii.pcently, Coel [ 16 3 obtained the critical speeds of 
uniform, isotropic undamped shafts using transfer FEH. The 
effects of shear deformation, rotatory inertia and gyroscopic 
moments were included in his analysis. It also covers 
multi-span shafts. In his work, two second order governing 
differential equations ^solved using linear polynomial , 
satisfying the continuity and completeness requirement led to 
ill- conditioning . To avoid this, these two equations were 
combined into a fourth order differential equation and 
cubic polynomial used. Although it was not possible to 
satisfy all the geometric and natural boundary conditions 
exactly in this approach, the results were quite 
satisfactory . 

1 .3 AIM AND SCOPE OF PRESENT WORK ; 

The present work aims at obtaining the critical 
speeds of flexible rotating shafts and studying the effect 
of various parameters on critical speeds . The relevant 
governing equations were derived and ualerkin FEM has been 
used to solve these more familier differential equations . 

The effects of gyroscopic moments, rotatory inertia 
and shear deformation are incorporated. A stud}? of critical 
speeds, neglecting the shear deformation is also done. Two 
second order governing equation approach is used in the presei 
work, except for the cases studied neglecting the exfect of 


shear deformation. 
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For the vibr 'tion analysis of shafts, carrying disks 
havin^ m sc larger than shaft mass, the mass of shaft is 
generally neglected. An attempt has been made to estimate 
the error committed due to this assumption. For this, critical 
speeds for a range of ratios of shaft mass to disk mass 
are obtained. Further, the effect of flexibility of disk is 
also studied. 

The effect of dissimilar shaft stiffness, caused by 
the unequal area moments of inertia of shaft section in the 
two principal planes is studied. The effect of dissimilarity 
in support conditions in the two planes is also included. 

A study of instabilities caused by internal damping 
and t he effect of bearing stiffness and damping properties 

on the stability aire studied. Hysteretic , os well as viscous 
forms of damping are considered. The combined model 
incorporating both the forms, is formulated based on the work 
of kelson [9 ]. The model is studied independently for both 
the forms of damping. Variation in instability threshold, 
with bearing stiffness and bearing clomping is also studied. 


CHAPTER- II 


THEORY A! ID FORMUTATT ON 


This chapter gives the derivation of governing 
differential equations and corresponding FSM formulation 
of the various cases of shafts studied in this work. 


2.1 ISOTROPIC SHAFT WITH NO DAMPING : 

Governing differential equations for transverse 
vibrations of rotating shaft, incorporating the effects of 
shear deformation, gyroscopic moments and rotatory inertia 
are derived from Hamilton's principle. 


Consider a differential element of length ds, as shov/n 
in Fig, (2.1). The kinetic energy of the shaft becomes f 8 j 

V4" (F+y 2 )d 5+ s 1 i D (et e 2 ) ds + 4 U o^s-/ i q 0 be ds 


o £ 

Here, m is the mass per unit length ; x and y are the 
translations along OX and OY axes 


; I D is the diametral mass 


moment of inertia per unit length ; Q x and are the bending 
slopes in YS and XS planes; 1^ is the polar mass moment of 
inertia per unit length and q q is the spin speed, assumed 
constant. Dot represents the derivative with respect to time t. 
In the above Eqn. (2.1), ’the first term is the energy due to 
translation, second term is the energy due to rotatory inertia, 
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the third term is the energy due to spinning of the shaft and 
the lost term gives the energy due to gyroscopic moments, 
derivation of kinetic energy is given in Appendix I. 

The strain energy of the shaft is 

1 1 
U = *jr / nl (0^ + © I ')ds + i / K GA (P 2 *- ^v ^ 
o v y o y 

.... ( 2 . 2 ) 

where El is the flexural rigidity ; K is transverse shear 
form factor, a constant depending upon the shape of cross 
section ; G is the shear modulus ; A is the area of cross section 
f x and j? are the angles of distortion due to shear deformation 
about X and Y axes respectively. Strokes represent the 
derivatives with respect to S. First term in Eqn. (2.2) 
represents strain energy due to bending and the second term , 
strain energy due to shear def orm-: tion. 


It is seen from the Fig. (2,1 Vthat total slopes are 


x = V'V 

y 1 = e v +/3 v 


.... (2.3) 


Using Eqns , (2.3), expression for strain energy U given by 
Eqn. (2.2) becomes 

U= ; i EI(Q 2 +e' 2 )ds + J \ t GA[ (x 1 -© ) 2 +(y» -© x ) 2 ] ds 
o ^ x y o 

(2. A) 

Thus integral to be minimised is 

I - CVJ) « 

t 1 



12 


*2 1 


0rl= i i 0 Q 

- x p °o ®xV 2 ei k 2+ ®; 2 ) 


- - K 1 GA [(x' -e y ) 2 +(y' -9 X ) 2 ]> ds at 


By Hamilton's principle 

t 9 

61 = 6 J (T-U)dt = o 


... (2.5) 


.... ( 2 . 6 ) 


Applying operator 6 to Eqn. (2.5) is quite straight 
forward, but fourth term needs little special attention, 


t 2 1 

6 q o 5 c ° 6 y ds dt 

1 t~ 


and is 


t ? 1 


' ' x p °0 6 <» x ) e y dt ds + / , I n e Jo as dt 


O t 


S I fi 6© © 

o P 0 x w y 


t 1 o p 0 x y 


1 t. 


ds ” Q ; / ^ D 0 6e A dt ds 

t 1 0 ■ * 


t 2 1 


2 1 

+ i l 5 °0 ®x 6e y ds « 


( o h °0 (-8 y 6e x + A 6e y) ds dt 

.... (2.7) 

Thus , applying the 6 operator to all terms in Eqn, (2,5) 
and rearranging, one gets 
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-to 1 . 

/" / ^[-m x+(K’GA(x< -© )) ] 5x 
t 1 o y 


+ [-I D VI D o e + (EISp' +K'G«(y'- 9 X ) ]6Q 3 


+ [-my + ( K' GA (y‘ -8.,))' ]6y 


+ [- X D Sy'S C o V< EI By)' +K ' a ‘( xL e y) 3 66 y 


ei e x ae x 


1 , 11 

-EI O' 6©. 


-k'ga(x / -e y )6 X 


l 


o 


-K'GA (y'-0 x ) 6y 


= 0 


The displacements 6x , 6©^, 5y and 6© y are arbitrary, 
(K f Ga (x* -0 ) ) - mx = 0 .... 

(EI 0y / + K'GA (x«-0 y ) - I D 9 y - I p D 0 Q X - C 


(K 1 GA (y' - © )) -'my = 0 


(EI 0 X )' +K , GA(y l - © V )-I D © X +I 


1 ) — o-r\ * jl 0 Q — 0 

x' D x p o y 

,1 


EI ©' 6© 1 = 0 ; K 1 GA (x* - 0 ) 6x j = 

y i o ^ *o 


o 


EI © x 6© x 


= 0 ; K' GA (y* - © x )’ 6y 


= 0 


o 


] ds dt 


( 2 . 8 ) 

therefore 

( 2 . 9 ) 


(2.10) 


( 2 . 11 ) 


( 2 . 12 ) 
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Equations ( 2 . V ) are the governing differential equations in 
IS plane, Eqn. (2.10) in YS plane and Eqns . (2.11) and (2.12) 

are the corresponding boundary conditions. These equations can 
also' be derived by Newton’s method [4 ] , (see Appendix II). 

Multiplying Eqns. (2.10) by i = 1, combining it v/ith 

Eqns , (2.9) and introducing the expressions 

z = x+iy and 0 = © + iQ v .... (2.13) 

y 

complex differential equations for isotropic shaft are obtained as 

(k‘ga(z' - 0 ))' - mz = 0 .... (2.14) 

(El 0’ / + K 1 GA (z' - 0 ) -I D 'i + il Q q ji = 0 

and for isotropic boundaries the complex boundary conditions 
become 

El 0 1 = 0 or 0 = 0 at x = 0 and 1 

(2.15) 

K 1 GA(z‘ - 0 ) = 0 or,z=0 at x=0 and 1 


The two equations given above can be combined into a single 
fourth order differential equation as 


El 




I ri (1 + ~ t ~) z il+ i I n C n z^+mz + 


m I- 


L B 


K‘ G 


p “o 


K* GA 


D T 


ml 0 
i — —2 — 2 z 


K 4 GA 


0 


(2.16) 

which Is same as that given by [ 4 ]. 

If shear deformations are negligible, the above 


equation becomes ; 
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El 


,iv 


I n 'z + i I G z" + mz = 0 .... (2.17) 


I) ~ *p o" 

The second and third terms represent the rotatory inertia and 
gyroscopic actions respectively. The above equation is same 
as that given in[ ]. 


2.1 .1 Hon dimensionalisation : 

The governing differential equations (2.1 A) and 
(2.17) are non-dimensionalised using the following parameters 


2 

1 


1 _ . O _ ® . rp 

, O - 1 


t 

T 


PT i c = Q 0 T 


* 


*n 


n T * 

, t - x . T - 

3”^ f •'•■vi T f 

0 0 


PAo 1 


4 


.... (2.18) 


El, 


Where end I are the area and area moment of inertia in 
o o 

bending at S = 0. 

Using these parameters, Eqns . (2.14) and (2.15) 

become 


(ii n (Z* - 0 )) - A n Me Re^ Z = 0 


n 


... (2.19) 


(j )' + 

n Me Re^ 


1 


(Z l - 0)- In Re 2 0 +i 2I n Re* Q 0 = 0 


Me Re 2 I n 0' = 0 or 0=0 
A n (z'~ 0 ) - 0 or Z = 0 

T 

E J n 2 o 
where Me = and Re - ^ 

K’ G A o l 

It may be noted that now prime (0 indicates 
differentiation with respect to S and dot (.) with respect to T. 


( 2 . 20 ) 


( 2 . 21 ) 


For uniform shaft A and I are unity , and 

n n J ’ 

Fqns , (2.19) and (2.20) become 

(s' - 0 i - Me Re 2 Z = 0 

( 2 . 22 ) 

0 + — ■* — k~ (Z f - 0)- Re 2 0 + i 2Re 2 Q 0 = 0 
Me Re ^ 


Me Re 2 0 7 = 0 or 0 = 0 
(Z 1 - 0) = 0 or Z = 0 

Similarly Eqn. (2.17) becomes 

Z IV - Re 2 (z" - 2i Q Z." ) + Z = 0 


(2.23) 


* • * • 


(2.24) 


2.1,2. PE?' I Equations : 

Ecu tions (2.19) are partial differential equations. 
Their solution is obtained by the finite element methoo . This 
method leads to simultaneous ordinary differential equations 
which can bo solved in a routine way. 


Descretizing the domain (0,1) by number of finite 
elements of length h (Refer Fig. 2.2) the variation of 0 and Z 

over a typical element can be approximated as 

/ N 7 (ne) 

z = a +bZ+ = [It \{Z) (2.25) 

, , ‘a (ne) 

0 e ^ = a'+b' 0+ = [IT J {_0 } 

where [N 4 * J and [ J are the shape functions ; {Z} 

^0 j(ne) are ^. he no( j_ a q parameters. 


and 





( 1 ) ( 2 ) ( 3 ) 

O' " — — - "O " O *■» 

12 3 4 

S 

s=o 


(c) 

j 




_h_ 

(e) 


H 

—o 

j 


— o . 

n (n+1) 


5 = 1 


Fig. 2-2 Finite element discretization 
and typical element 


1c 


substituting Eqns , (2.25) in Eqns. 
obtained as 


(2.19), residues are 


rZ = < A n ( z(e) ' - 0 (<O ))' - A n Me Re 2 Z (e) 

= U n «f (e) ’ )+ — — , ,. n (Z {e )' - 0 (e) ) .... (2.26) 

Me Re^ n 

I n Re 2 0^ + i 2I n Re 2 0 0^ 


Minimising these residues by the Galerkin method, 
one gets after integration by parts, 
h 


1 h 

J Ijf A n (Z (e) - 0 (e) ) dh +Mc Re 2 ' A n N 2 Z (e) dh 


= N 2 A n (Z^ e)/ - 0^ ) 


h 


o 


4 lj n ^ (S)/ db — o { N? A^(Z (e)/ - 0 (e) ) dh 


j N 


0 


h 


Me Re 


O ^ 1 V • Jrx \ 

2 0 1 n 


+ Re 


2 “ M 0 


h 

.2 r t \ T 0 


S N? I 0^ dh - i2 QRe 2 f. N? I dh 

M •*» ** U i ** 

■h 


N 


$ I 0^ e ^' 


i n 


.... (2.27) 

hatching the above equations, it is seen that compatibility of 
Z and 0 at nodes and completeness of Z, z! , 0 and 0 f is 
required. A linear polynomial i.e. C Q element for Z and 0 
satisfies the above requirements, but it leads to ill -co ndit ioning 
[16 ] . So a cubic polynomial i.e. element is chosen for 
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Z and 0 both, v.ith this, fewer elements need be used for 
same accuracy . Since N 2 = N^, the superscripts are dropped 
now onwards . 

Substituting for and 0^ e ^ and arranging Eqns . (2.27) 

in matrix form, one gets 

h , (ne) h (ne) 

j / ‘ n (NJ l> T Jdh (z? - ;a {n*> [Njdh^l 

0 0 


+ Me 

Re 2 

h 

; 



O' 

Me 

Re 2 

*h 

! s 



-0 

h 



+ / 
0 

‘"'n 

im 

- 2i 

Me 

4 

Re 


.. (ne) 


n 


n 


= {NJ A n (Z^ e)f - 0 (e) ) 


(ne) 


h 

o 

(ne) 


anc 


n 


(ne) 


4 


h 


(ne) 


h 


. (ne) 


- 2i Me Re / I n {N] [N J dh = {NJ Me Re* In 0 


( 2 . 28 ) 


These equation can be arranged in, the form 


McRe 2 [MA]^ e ) [0 ] 


w — , 

{zi A e ) 

Q 

-i 

To ] [o ] 


— —1 

IZ} 

[0] McRe 4 [MI]^ e ' > 


iiJi (ne) 


[0 ] 2I-icRe^[MI ] ( e ^ 


$ 0 } 


"[KA]<«) -C Ia, ] (e) 


_ {Zi (ne)' 

» 

FBS} 

-[KU] ( e ^ T McRe 2 [KA ] ^ e ^ + [MA ] ( e ^ 

•— — 


)0j (ne) 


f FBM } (■ ne ^ 


or [M] (e) ?qi (ne) - 10 [G] (e ^ < 4i (ne) + [K ] (e) *qj (ne) - }FB$ (ne) 


(2.29) 


where 


[ KA](°> 
[KU]( e > 
r 

[KI] 

[MA]^ 

[ ia] (e > 


h 

S 

o 

h 

J 

o 

h 

; 

o 

h 

S 

o 

h 

J 

o 


n 


n 


n 


A 


n 


n 


{N 1 } [N' J dh 

{N 1 } [N J dh 

in'} [ n'J dh 

{N} [Nj dh 
{.Nj [Nj dh 



{ FBI i \ 


(ne) f-McRe 2 In $ 


0 


McRe 2 I n 




0 



.... (2.30) 


— ( & ^ 
it may he noted that here all matrices [ M] v 

[G]^ e ^ and [K are symmetric. 

Equation {2* Zb) Is the governing equation v/hen shear 
deformation is negligible, -tyiis is a single fourth order 
differential equation. Its PEtl equations are obtained in the 
routine way and are 
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i’qj (ne) -io[G ]< e > ^ qj< ne >-, [K^ e ) { q } (no >= { F£.f ne > 

(2.31) 

where , 

(e) h ^ 

[M ] = ; {n} |_ NJ dh + Re 2 j fw'J I N'Jdh 

o o 


is the 


is the 


elemental mass matrix, 

r ,( e ) o h 

[G ] = 2Re 2 /jKl'U N'Jdh 

elemental gyroscopic matrix, 
(e) h 

[K] = /[N"}lN"Jdl, 


is the elemental stiffness matrix, 

(lie) (ne) 

is the elemental nodal displacement vector 
and 


ih 


(ne) (e)« 

{FB^ = -{H} 2 


+(N'j Z^' 


h 


-[n] Re 


2 pi 


-[N}2Re 2 0 # 


is the elemental nocial force vector. 


It may be noted that here element is needed, and 
all matrices are symmetric. 

2.2 DUAL SHAFT • 

Non-circular shafts or circular shafts with key ways 
and slots will have different stiffnesses m different 
directions, stiffness in any direction can be determined , once 
stiffnesses in two principal directions are known. Thus if 
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am Ig are the principal moments of inertia m principal 
directions 1 and II, then moments of inertia in X and Y directior 
aro, Fig. (2.3 a) , [17 ] . 


I 1 +I 2> 


I 1~ I 2n 


xx 


yy 


xy 


( — _) + ') cos 2¥ 

L,+In I^Io 

_ (-~-^) cos 2¥ 

i. -i„ 

( -' - g- — ) sm 2^ 


.... (2.32) 


Such shafts are referred as dual shafts m literature. 

For such shafts the moment curvature relations are [18] 

- - t EI xx 9 x + EI xy e ; ] 

.... (2.33) 

y L yy y *y x J 


Shaft is rotating with angular velocity q . X and Y 
are the fixed axes i.e. they do not change, «t any time the 
principal axes make an angle 0 Q t withX,Y-axes Fig. 2.3(b), 


Using the Eqns . (2.32) (¥=- fi Q t) and (2,33) and 
following the procedure of Appendix II, the equations of 
motion for uniform dual shaft become,, 

For bending in XS plane : 


VJ GA (x> - 0 ) - mx » 0 


E( - cos 2 ° 0 t)G(J- sin 2 0 Q t) 0. 


Vi 

2 


I- -I, 


.... (2.34) 


+ k'gA(x'- ©..)- In eL - 0 = 


0 



Fig*. 2.3 (a) Position of frame(XYS) relative to (I IIS) 
(b) Position of frame (IHS) relative to (XYS) 



?J\ 


For bending m YS plane : 

K 1 Gj* (y f - © )' - my = 0 


I„+I, 


I* -Ip i -i 

E( '~"T" + — cos 2 0 o t ) Q x" E ^ V ' " sin 2 V 3 ) 6 y 

+ K ' GA (y'- 9 x) - I D ®x + !p "o Sy = 0 

Multiplying Eqn, (2*35) with i = J~-1 and combining with 

Eon. (2*34), using the notations z = x+iy ; © = 0 + i© and 


(2.35) 


© = © y - i© x , the equations of motion become 
K 1 GA ( z 1 - 9) 1 ~ m’z = 0 

I 1 +I 2 ./ Ii-I? „ 

E( - lg- ^ -) 9 - E cos 2 D 0 t 5 

I -I 

-iE ( — -g— -) sm 2 C Q t 5^ + k'ga (z‘ -9) 


- ! D 9 + Up 2 0 © = 0 


.**. (2.36) 


Combining second and third terms of second equation, one gets 
K* GA (z 1 - 9)‘ - mi = 0 


E(^l^)© 1 '- E( —g— ) S'e l2 £ o t + K*Ga (z* - ©) 

- I D 9 + H p D 0 9 - 0 

-i 0 t 

Multiplying Eqn, (2.37) by e ° 

k'ga (z l - 9)' e" 1 V 1 - mz e _i C o t = 0 


..*. (2*37) 


E (il^ 2 )Q"e 1 V e 1 ^ +k'gA(z'- ©)e~ iQ o^ 


• * * » ( 2 « 38) 


-Ir, 9 e” 1 Q o t + il„ Q, © e” 1 Q 0 t 


0 


Equations (2.36) have coefficients which ere functions of 
time. Tnis time dependence can be eliminated by defining 


£ = z e 


-i B o t 


-i fl t 
0 = 9e 0 


0 = 5 


Q t 
o 


(2.39) 


Thus Eqns . (2.38) become 

K l CA (V -0)‘ - mV - i2m t- * 0 

I, “Ip T if I y 

E(-V^)^ - E(-l^) 0 + K'GA (S' -0) 

.... ( 2 . 40 ) 

_I D y - 12 I d Q 0 8i + I D 0* t 

* “p O o<* “ *p B o * = 0 

Taking I D = ^ I , as done in [4 ], the above equations become 

i O 

iV GA ( S* -0) - mV -i2m C 0 S+m 0 q S = 0 

E (b^2) f(". £(3^) tfV'GAU'-ltf) .... (2.4D 

- Id * - Id 4 81 ■ 0 

These are differential equations in complex variables 

v l7 V ^ ~ e v +lQ '% * ^ ~ ®rT iG *S “•* ( 2 *^ 2 ) 

w iere ^ n are the deflections along I and XI respectively 
and 9^, 9^ are the bending slopes m IS plane and IIS 
plane -respectively , 


3epr,rvitin_ real and imaginary parts, one gets 


to 


e n )'H n 0+ 2mB o n o+ m 0^ s 0 = 0 

EX 2 e" + K , GA( <i ' o -9„)-i D e n -i D B z 0 e n = o 


K'G/,(ii 0 - )' -m r ’ 0 -2m o Q % Q +m c^ n 0 = 0 

El, +K'cA(r/- Sj) -I D «ij -I D s^ = 0 


(2. Vl) 


Equations (2.43) correspond to IS plane and 
Eqns . (2,44) correspond to IIS plane. 

It may be noted that Eqns, (2,43) and (2.44) differ from those 
given by Dimentberg [ f ] . It is because of the definitions 
of 1^ and Irpfwbfch differ from Dimentberg' s definition. 

By observing Eqns. (2,39) > it may also be noted that, 
the natural frequencies obtained using these conversions, 
will be the natural frequencies (\ ; ) in I II S frame , 

To obtain the natural frequencies (4) in XYS frame, the 
following relation is to be used, 

4 — 4 + q q , « « , ( 2 , 45 ) 

Two second order equations of Eqn. (2,43) can be 
combined into a fourth order equation and same can be done 
for Eqn. (2,44), These two fourth order equations are same as 




tliO'-.c given by Dimentberg [4 ] * which were obtained m a direct 
\h’Y without deriving the second order equations. These second 
order equations are of great advantage m FEU formulation, 
as all boundary conditions (geometrical and natural) can be 
incorporated exactly. When shear deformations are negligible, 
these tvo fourth order equations become for bending in 
IS plane, 


El 


2 d 0 v -i D c^ 


o 


r )+m( t-2 CV ^ O = 0 

o' v ^o o o o u o / 

.... ( 2 . 46) 


and for bending in IIS plane 


EI 1 + °o"o> *»> d 0 +2 




(2.47) 


2,2.1 Non-dimensionalisation : 


The non-dimensional parameters used are 
*>o n 




T~ ’ r> ~ 1 


; S* f ; T 4i 
1 T 


Q = DqT . T = ; C= jL , C d = ^ 


, , , , (2 ,48) 


where k is the radius of gyration for mass moment of inertia 
and is the radius of gyration for area moment of inertia , 
With those the non-dimensional form of Eons. (2.43) and 


(2.44) are 


2S 


J ~ J 


( Zy- ©„)- Me Re^ ( t, -2 QT)~ ) = 0 


McRe 


(2.49) 


Q ‘t) + Ctl - e n) " c d ^ Re^(©^ C ^ ©„) = 0 


and 


(V- ©£ /- Me Re 2 (n + 2 0 -■E.Q 2 ^) = 0 

lie Re 2 ©" + (n* - )-C d Me Re^ (e + o 2 e^ 


.... (2.50) 
= 0 


2 ^"1 

where Re,. = — *— , is the rotatory inertia co-efficient and, 
A1 

g 

Me = — r— is the material constant. 

KG 

Similarly Eqns . (2.46) and (2.47) become 


E, IV ~ C C d Re 2 (*." + 0 V)+C (i -2 0"- D 2 E, ) = 0 (2.51) 

and 

rj IV - C d Re 2 (t/+ c 2 Ti") + (ri + 2 ®fc,~ Q 2 ^) =0 .... (2.52) 


2.2.2 FEM Equations : 


Proceeding as in section (2.1 .2) finite element 

(ne) 


equations for dual shaft are 

(e) (ne) (e) . (ne) (e) (ne) 

[K ] i q‘3 +[&] , +U3 {q] 


=5.fb} 


• # « • 


(2.53) 
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where 


[il (e} ] = 


t G^] = 


He Reij; M^] 

[0 ] 

[0] 

[0 ] 

[o ] c 

^Mc Re^ 

[Hi] [0] 

[0 ] 

[o ] 

ro ] 

Me Re 2 [ ] 

[0 ] 

_[o ] 

[o ] 

[0 ] 

C d Mc Re 

elemental mass 

matrix 

f 


— «* 



— 

[0 ] 

[0 ] 

-2 0HcRe^[ M 1 ] [ 0 ] 

[0 ] 

[0 ] 

[0 ] 

[o] 

2 QMcRe 2 [M^ ] 

[° ] 

[0 ] 

[0] 

A 0 ] 

[0 ] 

[0 ] 

[0]_ 

elemental gyroscopic i 

natnx, 


C[ K 1 ]-McRe 2 

olMp) 

[K 2 ] 

[0 ] 


r >< (e) ] 


McReV 

(-7T 1 M 


[0 ] 


[0 ] 
[0 ] 


+(1+C d McRe^ Q 2 ) 


e>ym , 


[i^]) 


K 1 -Me Re 2 q 2 


[K 9 ] 


[Mil 


( McRe 2 [K^l 
h 


L 


+(1+C d McRe^ 
• 2 2 ) [%]) 


is the elemental stiffness matrix, 


tu (ne) 

(ne) - W ne) 

{rj] ( ne ) 

l&i (ne) . 


is the nodal displacement vector of the element, 

{N} - 4 e) ) 




(ne). 


’1 C>n 


McRe 

-IT- 1 l N ] G n 

{%} (r> (e) - (J e) ) 


McRe* {NJ 

V» 


T) ' 'o 

h 


o 

h 

o 

h 


is the nodal force vector of the element, 
and 

h 

[I-m] = S m [Njdh 
[^] - ; IN'] [N'jdh 


.... (2.54) 


h 


[K ? ] = - / tti 1 } L N J dh 

c ' o 

f V ^ 0 ^ / \ 

Here [ M ] ^ e ^ and [k] are symmetric matrices, and [G 
skew -symmetric matrix. 

Similarly FEM equations for a dual shaft without 
shear 'Sqns. (2.51) and (2.52)„become 


is 



where, 


r M 


] (e > ^ + 2 [c- ] (®) { c; 


'C([ Ml ] +C d Re 1 2 [K 1 ]) 

r K ] ^ e ^ : 

[ o] 


is the elemental mass matrix, 


(ne) 

[F3] (2.55) 


[ 0 ] 

( [ i'L,] + C d °e^ [ K^j] ) 


[G] 


(e) 


[ 0 ] 

2 [ Mj 


-2C [Mj 

[0] 


is the elemental gyroscopic matrix , 


[K] 


(e) 


[K 2 ] +C 0 2 (C d Re 2 [ly] ■[&,] ) 

[ 0 ] 


0 


[K 2 ]+ Q 2 (C d Re 2 [IL,] 
'[Ml] ) 


is the elemental stiffness matrix. 


[T 3 ] 


( ne) 


-l v ( e ) 


■m 


jh 


?N] ri 


(e)« 


o 

h i y 

+ {K'jrr e - ) 
o 


n 

o 

i/ih 


h 


+CC d Re^N}^ ( e ^ +CC d Re^ rf $11$ ^ e ^ 


+C d Re 2 [N?ri 


(e) 


+C d Re 2 c 2 {N^( e ) / 1 

o 


0 


is the elemental nodal force vector, 



is the elemental nodal displacement vector 


snd 

h 

[;l] = i inj dh 

O 

[ll,] = / f N* J [ N 'j d h (2.56) 

o 

[Kp] = ; fN fl } L K l, Jdh 

o 

For both the cases element has been used as 
discussed in section (2.1 ). Here [K and [h ^ are 

symmetric matrices, and [g is unsyrometric matrix. 

2.5 IoPTROPIC SHAFT V/ITH DAHF'ING : 

The Triction offers resistance to motion and this 
usually leads to damping of natural vibrations. The friction 
arising between rotating and stationery parts ,such as bearings 
is c- lied as external friction and damping arising due to 
•this is called external damping. The friction arising within 
the rotating shaft due to resistance of the particles of 
material is called internal friction and the damping due to 
this is called as internal damping . 

External damping can be incorporated by modelling the 
bearing as dash-pot model with a viscous damping - o-oeificient C 




ana e'gr.tin^ the shear force to t,he dampin^ force, while 
npj 'lying the boundary conditions , 

To take internal damping into account, two forms of 
damping moo els are taken. In the first, it is linear 
velocity dependent viscous model, for which the constitutive 
relation relating the axial stress c to the axial strain 
e s is [ 9 ] . 

°s = \ E ....(2.57) 

where 71 ^ is viscous damping coefficient. 

The other model is frequency independent hysteretic 
dc- mping which is expressed by the loss, factor h which is 
related to loss angle 7^ [19] by 


tan 7 n = 77 h 

The constitutive relation for this case is [ 9 ] 


. (2.58) 


- E [ (-=hr) e s + < 


n 


H 


JT+n 


\JvPn. 


■K 1 


. (2.59) 


H 


Combining the two models, Eqn. (2.57) and Sqn. (2.59) 
one gets 


71 . 


=E [(■ 


JiTfi. 


.)e + (71 + 


H 


Juri 


-)e ] 
‘ s J 


.... ( 2 . 60 ) 


H 



Koferrin^ to Fig. (2.4 ), z (s,t) is the distance of 
neutral axis of cross section, X is she speed of precession 
C Q is the spin speed. The normal strain in the differential 
area an , at a radius r , is 


h - - -rcos( O -Mt0 


(2.61) 


where o is the perpendicular distance of dA from neutral axi 
From Fit,. (2.4), moments about x and Y axes are 

2n r 


M = S /°[ y + r sin Dt]o g dA 
x oo 


.... ( 2 . 62 ) 


2n r 

K = - j / [x+ r cos D t J c s dA 
y oo 


Using Eons. (2.60) and (2.61), moments from 
.on, (2.62) become *[ 9 ] » 


I r v 

A 


ri-r 
= £j1 


M 

y 


-(- 


1+n 


H 


/1+n. 


-) 


( Q!hj+ 


n 


H 


■) 


n 1 +tu 

(^=r> ( =f 

i 


) 


1+n 




x 


J 


+ El 


1 

< 

o 

_ l 


' yH 



* «( 

0 


X 



_ — 


.... ( 2 . 6-0 


A l'> 


A*' 


2,'H 


Using the above Eqn. (2.63) in Eqns . (A-*) and (#=M of 
Appendix II, 'the governing differential equations are 


obtained as ; 


1 + 71 , . 

El [ U r iv + ( . 


'H 


+ n v D J y 1V +T V ^ 1V ] - 


1+71 


H 


F 


V 


D 


+ 71 . 


H 


, // 

I D v + mx = 0 
p o J 


(2.64) 


El [ 


for bending i n XS plane and 

1 i-Tl,. 71 . 

r 1 - y 1V - (—£- + b v Qq ) x lv +T7 y y lv ] 


.. fi 


H 


-I D y 


H 


+ I c x"+ my = 0 
p o J 


for bending in YS plane. 


2.3.1 Non -dimensional is at ion : 

The non-dimensional form of Eq ns . (2.64), are obtained 
following the procedure of section (2.1.1). The non- 
dimensional parameters used are, 


y ~ . y ~ X . t = — f r ■ 

x - T ' Y “ 1 » 1 "7 ' 


Q = 


* 


o 


71 


7 *, 4 


= ^ ; 


V . m* _ [ PA1 
’ 1 'v El 


.... (2.65) 


k 1 

Tr XX p, 

K xx = El ’ XX 


c xx 1 
El y 


Then, Eqns. (2.64) in the non^iimensional form, are 


(-XX) X IV + (XX. + 71 o ) Y IV + h VN 

V^h 2 JH ■ 



Re 2 (i l> + 


2 0f) 


+ X = 0 


; / 


1+rj u 
( ul) Y 


where Re' 


IV 


C— — + r \„- Q ) :: IV + T ‘ 


fur,- 


V 




. (2.66) 


O •> // . // * * 

- Re^ ( Y -2 C X j + Y = 0 


I 

-H.L 


2.3.2 FEi. Equations : 

Following a procedure similar to the one explained 
in section (2.1.2) in detail, FEM equations are obtained and 
expressed in matrix form as f ollov/s . 

/ \ f \ ( e ) (ne) (e) (ne) (ne) 

[M] (0) 13] (ne + [ °] Vi +[K] w 


(2.67) 


where 


[M ] (e): 


[Ml] + Re [ K ll 
[0 ] 


[0 ] 

[M-,] +Re 2 [r.J 


is the elemental mass matrix, 

n [K 2 ] 2 ORe 2 [K^ 


[C ] 


( S ). 


-2 £Re" [It, ] 
is the elemental damping matrix, 




[K] 


(e) 


r 1+71 H r 1 

( -) [k 2 ] 

(— 

, p 

/t+ng 


-(-V- +*V M o )[ k 2 ] 

1+T) H 
^ ri 


J 


+n VK C - 1 [*7 


) [K 2 ] 



is the elements! stiffness matrix, 





.fxj 

(ne) 

(ne) 



(ne) 


jYi 



is the nodal displacement vector of the element, 

(ne) 


{FB} 


(ne) 


* F d 

j F z) 


(ne) 


is the nodal force vector of the element, 
where 




(ne) = |- _ |>j ij ySef 


h 


1 


i+n 


H 


n 


H 


+ (— '"VU 




+ ”v 13 C ) [ - fK} Y^ e)/ " M N 'h 


h 


-(e)" 


h 


+ ^VN l 


r- X^ 


/// 


■h 




+ Re^ [ ^ X 


Me) 


i ih 


{N 1 } 


{%} 


h 


y(e) 


o 

I ,h 
0 


{ F ?1 


(ne) 1+r , H 


(±U1 ) r - {N} Y 

fun 


(e) 


ut 


•h 


+ [N ; ] y^ U 


h 


o 


H 


rj„ -fpM/i 

- (-S- * r ' W > [ - l N j ^ ' 

^ I 

•h ih 


h h y ( e ) 


t/ ih 




o 


+n 


VN L 


W 


n y( e ) /r/ + ^n 7 } 


, w h . (e ) 1 

Re2[*N\Y (e >' - 


o 


o 

h 

o 



and 


[Ml] = 


[K-] = 


[Ko] 


h 

I {N] [ N J dh 
o 

h / 

j {N } L N'jdh 
o 

h „ 

{N }[ N Jdh 
o 


( 2 . 68 ) 


[M] (e) is symmetric ana matrices [C] 1 ' 6 ' 1 and [K] v J are 


Here also element has been used Matrix 

-i(e) rvi( e ) 


skew -symmetric , 


CH/’.FTER- III 


RESULTS aKD DISCUSSION 

\ 

Computer programs in FQRTRAI’ IV have been developed 
to determine the natural frequencies and critical speeds, 
for all trie cases of shafts, formulated in Chapter 2 . 

3.1 METHOD OF SOLUTION : 

^ * 4. • r t t l( e ) r r l(e) r r i(e)and f [K]^ e) , 

Element matrices, [Mf, [G J , L C J 

were calculated using Gauss quadrature with four Gauss 

points. Those matrices are assembled in the usual way and 

boundary conditions applied. 

3.1.1 Isotropic shaft without damping : 

Assuming a hormonic solution of the form 

n iAt ( 3 i ) 

{ 0 ] * e * * * * ' 

ond applying to the assembled form of non. (2.2S) the 
lollop In, eduction is obtained 

( n ) 

(E\ 2 [ K ] [G ] + [KDfq^ =0 .... (3.2) 

For critical speeds in forward precession, tne criterion is 

* • * • ( b .3 ) 

and for backward precession it is - 0 * 

Applying these critaria., ^qn . (3.2) becomes 


W 


o 


B 2 ( [M]T > ^ q L 


0 0 0 0 


(30) 
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wnore minus sign is to be used for forward precession and 
plus sign for backward precession. The eigenvalues D 
are the critical speeds . The above eigenvalue problem is 
solved using RIG library subroutine F02LEF. 

3.1. 2 Dual shaft without damping : 

The method of solution for determining the critical 
speeds of dual shaft is given in this section. Governing 
uifferentiol equations for this case are in reference frame 
I 113 . The natural frequencies (X 1 ) in this frame of reference 
ere related to natural frequencies (X) in frame of reference 
XYS , by Lqn. (2.45). It is to be noted here that in the 
bqn. (2.45) is the spin speed in dimensional form, which is 
later non dimens ionalised as q . 

Noting that, the critarion for critical speeas to occure 
in forward precession is X = c , in reference frame XYS , 
the critarion in I IIS frame becomes 

x' = 0 (3.5) 

Considering a dual shaft without shear, Lon. (2.55), the [ K] e 
matrix can bo expressed in the following form. 

(e) (e) 2 

[k] = [kd] -o [md] 


• • * • 


(3.6) 


k 2 


where 


[KD] 


(e) 



[ 0 ] 

[k 2 1 


and 


[ Ml. ] = 


C ( C M 1 1 - C d Re2 [ Kl ] ) [0 ] 


[ 0 ] 


-[^]“C d Re* [Zj 


Then ^qn, (2.55) can be re-written as 


{q^ ne ^ + [ g] ^4^ ne ) + ( [ KD] D 2 [md] )^q j ^ ne )= fFB^ ne ^ 


assuming a solution of the form 

I01 ln) -W D) ^ 


.... (3.7) 


(3.8) 


and applying to assembled form of Eqn.(3,7) and noting the 
condition Eqn. (3.5), the eigen value problem becomes 


[KD]^q o | (n) = Q 2 [MD ]^q^ (n) .... (3.9) 

in v/hich, the eigen values represent critical speeds in 
forward precession and Eqn. (3.9) is solved using II vG 
library subroutine F02A.EF. 

For critical speeds in bac Inward precession, the critarion 
in aYS frame of reference is 


X 


0 


* » t < 


(3.10) 


A3 


Then, from Eqn. (2 . Ad), the criterion in I Ho frame 
becomes 

X 1 = ~ 2 Q .... (3.11) 

assuming solution of the form given by Son. (3.8) and 
applying to the assembled form of Son . (3.7), the eigen value 
problem to be solved becomes 

[ KO] m = 0 2 (A [K ] + 12 [G ] + [MD ] ) 5 a 7 .... (3.12) 

O ' ' os 

where the eigen values represent the critical speeds in 
backward precession. Nov; as one matrix is complex, the above 
eigen value: problem is solved using the NAG library 
subroutine F02GJF . 

The critical speeds of a dual shaft with shear 
deformation effect are also obtained using a similar 
procedure . 

3.1,3 Isotropic shaft with damping : 

The nature 1 frequencies for this case are obtained 
solving Bon. (2.67), for a range of spin speeds, using the 
method of solution giveb below. 

The Eqn. (2.67), after assembling and applying 
bound: ry conditions becomes 

[ K] ^l (n) h C] 5 q] (n) + f K] = 


0 


wilier c;.n be re-written as [20 ] 


[ 0 ] 

M 


whore 


[M] 

[C] 

m = 



- Cm] 

[o] 

[o] 

[E] 





Assuming n solution of the form 


0 


• * * * 


(3.13) 


5 p ( = |P 0 j e xt .... (3.14) 

where X = X,, + iX^ is complex, Eqn . (3.13) becomes an 
eigen value problem of the form 

[0] [I] 

[K f 1 [ M ] -[K ]" 1 [C] 

The above ei p en value problem is solved using HAG library 
subroutine F02 aFF. The eigen value? appear as complex conjugate 
pairs, with positive or negative real parts. The positive 
real port indicates an instability. 

3.2 COIhCXiSMCE 1TUDY : 

To study the convergence , results (critical 
speeds) were obtained for a steel shaft of uniform circular 
section, using 4,6 and 6 elements, by solving Eqn. (3.4 ), for 
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hinged -hin_,ed, fixed-fixed and cantilever boundary conditions, 

T? 

The value ol’ material constant (he = — — ) for such £ steel 

K' G 

shaft is 2.9. The value of rotatory inertia co-efficient 
(Re - used is 0,05, which represents a shaft having a 

length, five times its diameter. 


The results are presented in Tabic-. 1, For hinged -hi nsec, 
boundary conditions the results are compared with the closed 
form solution given by Dimentberg [ 4 ] . For other boundary 
conditions, comparison is done with the results obtained by 
Gocl [16 ] using transfer FBM with a single fourth order 
governing differential equation. 


„ith an increase in number oi elements used, the 
results converge from above and the error is larger in higher 
modes, as is to be expected by FBM. The critics! speeds 
obtained using 8 elements -were found to give accurate results 
for all the boundary conditions, for the number of mooes studied 
Hence 8 elements were used for all the cases studied in 


this worh. 


TaMe 1 : Convergence otudv of 

Critical Speeds 




(ae = 0.05, Txc = 

O O \ 

<— • y J 




Hinged -hinged 0 

hi ft 


node 


Ho. of elements 


— 


_J ^ 

5 

8 

— nxact [_ 4 J 

1 

9.6387 

9.6383 

9.6383 

9.6383 

2 

33,8967 

35.8793 

35.8774 

35.8768 

3 

72.3132 

72.1620 

72.1429 

72.1370 

4 

113 .6066 

112.3852 

112.2958 

112.2660 

5 

136.3399 

153.5009 

153.2122 

153.1123 



Fixed -Fixed 

She ft 


Mode 

5 

No. of elements 


T ~ 

0 4 

6 

E 

5 s r i6i 

1 

13.3408 

19.3379 

19.3377 

21 .660 

2 

47.0607 

47 .0430 

47.0367 

54.682 

3 

81 .5682 

81 . 3770 

81 .3501 

94.827 

4 

120 .4545 

118.9859 

118.8833 

136.970 

5 

161 .4977 

158.0973 

157.8044 

- 



Cantilever 

Shaft 



l 

No. of elements 



i'KJUC? 

1 25 — 

6 

8 


1 

3.4770 

3.4770 

3.4770 

3.477 

2 

20.4524 

20.4492 

20 . 4490 

20.441 

3 

52.0531 

51 .9979 

51 .9915 

51 .902 

4 

90.4899 

90.2031 

90.1609 

89.801 

5 

132.9002 

130.9754 

130.8163 

- 




hi 


3 . 3 lA'JiLTIC SHAFT , ITVJUT p I1FIIIG : 

-'.3.1 Jffect of gyrosccoic action md rot: torn inertia : 

Gyroscopic action term is the [G ] matrix in 
■tqn. (2.29 j. neglecting this term, the ecu- lions become as 
those of Timoshenko beam. desalts ( critical speeds ) for 
forward precession for v rious shafts are compared with those 
of Timoshenko beams for Ke values of 0.0125 and 0.05 in 
T: ble 2. 


Results for hinged -hinged and cantilever shaft match 
well with the exact results of [4 ] and [ 15 ] , but there is 
large difference for fixed-fixed shaft. This can he explained 
by notin, th t Goers [16] results are based on single governing 
.fourth order differential equation, whereas present results are 
based on two second order differential equations. In Goel’s f )6 ] 
approach one can satisfy only one geometrical boundary condi- 
tion (GbU) Z = 0 exactly. Other geometrical boundary condition 
of slope and two natural boundary (NBC) conditions on shear 
and bending, moment can be satisfied only in an approximate way 
where: s in the present, approach all boundary conc.itions 
(geometrical and natural) can be satisfied exactly. In cese of 
hinged -hinged shaft ,<3 oelfe approach satisfies u" r o GBC exactly 
and two NBC approximately. And for fixed -fixed shaft it 
satisfies two GBC exactly and other two GBC approximately. 

It is seen that approximate application of two GBC in case of 
fixed— fixed shafts leads to a large eri Oj. . 



Table 2 


.ytudy oT ryroccop ic and rot-tory inertia 
effects on critical speeds : 

(He = 2.9) 

Hinced-hi nred 


1 lode 

j Re 

= 0701 23 


. l6— 0 , 

.07 


J -imosnenko i 
J beam ; 

) Shaft 



[ Shaft 

Willi 


mm 

wsmmm 

Present 

1 

9.84 

9.85 

9.86 

9.43 

5.64 

9.64 

2 

39.02 

39.26 

39.25 

33.80 

35.89 

35.88 

3 

Cb . 58 

87.77 

87.70 

63 .47 

72.22 

72.14 

4 

151 .14 

154 .81 

154.43 

103.11 

112.66 

112.30 

3 

231 .11 

- 

238.43 

141 .54 


153.21 


Fixed -fixed 


Mode 

i Re 

= 0 .01 25 


' Re 

= 0.03 


oTimoshenko 

& Shaft 



I beam 

mwsm 


Present 



Present 

1 

22.11 

23 . 34 

22.15 

19.21 

21 .66 

19.34 

2 

60.02 

61 .29 

60.42 

45.04 

54.66 

47.04 

3 

115.42 

119.39 

116 .97 

76.53 

94.83 

81 .35 

4. 

186.47 

195.84 

190.53 

110 .85 

136.97 

118.88 

5 

271 .53 

— 


279.92 

146.96 


157.80 


Cantilever 



! Re = 0.0123 

Re = 0.03 

Houe -j 

[Timochenko j 

! Shaft 


1 b;-:n ] 






3,51 

3.51 

3.51 

3.94 

3.4s 

3.48 

21 .82 

21 .93 

21 .93 

19.23 

20.44 

20 .45 

60.32 

61.03 

61 .01 

4 / , 1 4 

51 ,90 

51 .99 

116.10 

118.62 

118.46 

80.15 

89.80 

90.16 

187.83 

— 

193.57 

115.91 

- 

130.82 


2 
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ncrunso the natural frequency of the system and the 
effect is signif icent in higher modes. This is because, the 
moments arising due to the gyroscopic action, cry to 
straighten tne snaft, and thus account for an increase in 
effective stiffness, which increases the n turn! frequencies 
01 the system. The gyroscopic moments increase with increase in 
curvature, -s the curvature is larger in higher modes, the 
effect is of appreciable magnitude in higher modes . The 
effect is seen to increase with an increase in lie. A larger Re 
implies a t hie leer shaft (i.e. a smaller l/d r'tio, where d 
is the diameter of shaft), which also leads to larger 
gyroscopic moments, and hence -n increased effect. 

It is also seen that critical speeds decrease with 
increase in Re value, and this effect is significant in 
higher modes, as expected. 


3.3.2 affect of shear deformation : 

Critical speeds in forv/urd precession are obtained 
by solving Eqns . (2.31), and are given Table 3 for various 
shafts, for two values of Re. Results with shear are also given 
in this table . 

It is clear from Table 3, that neglecting shear 
deformation effect leads to an over -estimation of critical 
sneedc and the effect is of appreciable magnitude for thicker 
shafts (i.e. higher Re), which is as expected. For Re = 0.05 



Tel: le 3 : 


Study of sheer deformation effect 
on critical speeds . 


(Me =2.9) 


Hinged -hinged 


Mode 

11 Re = 0 

.01 23 

\ Re = 

0.03 

ithout shear 

I With sheer 

{ fit}} out 

shear fcWith shear 

1 

9.88 

9.86 

9.99 

9.64 

2 

39 .61 

39.25 

41 .60 

35.88 

3 

89.57 

87.70 

100.84 

72.14 

L\ 

160.53 

154.43 

203.81 

1 1 2 . 30 

5 

253.98 

238.44 

402.69 

153.21 

6 

518.80 

338.57 

1092.45 

193.89 

7 

739.01 

453,71 

- 

234.20 

8 

933.83 

589.84 

— 

275.21 


Fixed -fixed 

fi Re = 0.0123 X 

Re = 0.07 

{Without shear 

With shear { 

Without shear \ With shear 


1 

2 

3 

4 

5 

6 

7 


22.40 

22.13 

61 .94 

60.42 

122.14 

116.97 

203.86 

190.53 

305.9^ 

279.92 

439 .68 

382.80 

592.05 

500.15 

847.90 

539.98 


22.73 

19.38 

65 .60 

47.04 

139 .66 

81 .35 

266.20 

1 1 8 . 88 

521 .53 

157.80 

1926.30 

197. 16 

- 

236.59 


277.08 


8 
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Table 3 • • • . (contcl. ) Cantilever 


Mode 

| Re = 0. 

0125 1 

1 Re 



5 With shear i 

\ V/ithout shear 

T ~ l 


mnnmnn 

5 

5 


1 

3.52 

3.51 

3.54 

3.^8 

O 

L— 

22.09 

21 .93 

22.99 

20.45 

3 

62.11 

61 .01 

66.72 

51 .95 

4 

122.55 

118.46 

150.87 

90.16 

5 

204.72 

193.57 

303.32 

130.62 

6 

310.38 

285.20 

679.33 

171 .84 

7 

441 .62 

392.04 

- 

212.56 

8 

593.75 

511 .67 

- 

252.90 
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only six critical speeds are obtained in forward precession 
for all the three shafts. This confirms the work of 
Dimentberg [4] , who showed by close-form solution that the 

number of critical speeds (n) in forward precession is 
limited by (for hinged-hinged supports), 

n i nRe~ •••• (3.16) 

Observing the results for fixed -fixed and cantilever 
support conditions, it seems that above relation holds good 
for these shafts also. 

3.3.3 Mul ti -span shafts : 

Three shafts with two, three and four spans used 
by Goel[l6] are studied and are shown in Table 4. All 
supports are hinged. Number of elements used in each span are 
shown in brackets. Critical speeds are presented in Table 4. 
Also presented are the results of rotor given by Gael 
using transfer FEM. 

With increase in number of supports, the stiffness of 
the system increases significantly. Hence, a significant 
increase in critical speeds is to be expected, and it is clear 
from Table 4. 

Except the results for first mode for three and four 
span shafts, the present results are comparable with those of 
Goel [16]*» These two results for the first mode are clearly 
in error as they violate the fact that increase in number of 
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Table 4 : Cri tiosl speeds for multi -sp-^n rotors 

( Re = 0.05 , He = 2.'-') 

Rote : 1) The numbers on top are the ratio of span length/ 
total length 

2) The numbers between two supports are the number of 
elements considered for that '.span. 


Two -span rotor; 


_ 0-3 

& ^ s 


0-1 


(S) 


Mod e 

Rotor 



Present 

Goel [16 ] 

1 

23.98 

24.95 

2 

69 .60 

70.62 

3 

95.87 

104.22 

k 

128.92 

138.36 

Three. 

-span rotor : 

ft. , Q-4 ...ft-g-fL 
Jr? (4> Ql (2> 



■ Mod e 

Rotor 

Present 

Goel [ 16 ] 

1 

60 . 1 0 

23.88 

2 

89.88 

86.40 

3 

113.42 

111 .91 

4 

160.22 

142.26 


& 


0-1 
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Four -span rotor : 


o ag o-g, ojs o-as - 
A C3) Au>fl (21 (2)5, 


Mode 


Rotor 

Present 

Goel 

1 

100 .26 

31 .69 

9 

t- 

121 .12 

112.63 

3 

128.78 

133.25 

4 

180 .64 

162.73 
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supports, increases the stiffness significantly. Ilorever these 
types of inaccuracies are quite common in transfer FEM. 
Otherwise results differ, as Goel [16] used a fourth order 
differential equation and present work uses two second order 
differential equations . 

3.4 EFFECT OF SHAFT MASS AND DISK FLEXIBILITY : 

As said earlier, most of the work in rotor dynamics 
neglects the mass of the shaft and the flexibility of the disk. 
Effect of these assumptions is studied in this section. 

To study the effect of shaft mass, a steel shaft of 
1000 mm length and 40 mm diameter with hinged -hinged supports 
is chosen. This shaft has a rigid disk of 30 mm thickness and 
90 mm diameter placed at a distance of 250 mm from the left 
end. For shaft, rotatory and gyroscopic effects are taken into 
consideration, i.e. governing equations are (2.31). Shaft is 
divided into eight equal length elements. Disk is taken to be 
rigid, but has rotatory and gyroscopic effects. Thus for the 
disk ut node 3, one has [4 ] 



(3.17) 
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where ra, in the mass of dish, I, and I are the diametral and 
a ’ a p 

polar mass moment of inertia of disk, and is the spin 
speed of shaft. Equation (2.31) and Eqn. (3.17) were 

assembled together and solved for critical speeds for a range 

% 

of shaft mass. Range of shaft mass is obtained by successively 
decreasing the density of shaft and maintaining all other 
parameters (length, diameter and E) constant so as not to 
alter the stiffness of shaft. Critical speeds for first five 
modes as a function of ratio of shaft mass to disk mass are 
given in Table 5. 

The critical speeds for massless shaft are obtained 
by reducing the order of assembled stiffness matrix to that oj. 
mass matrix of disk i.e. 2x2, using static condensation [27], 
For static condensation, see Appendix III. For forward 
precession ,only one critical speed exists and is sho*n in 
Table 5. This is same as the close form solution given by 

Too [22]. 

It is seen for this shaft and disk, in the first mode 
there is an error of 156 */. for a shaft to disk mass ratio of 
6.564 and 30 */. f or a shaft to disk mass ratio of 0.844, 
Moreover massless shaft assumption gives only one moae m 
forward precession. Thus one should not neglect the shaft 
mass, as error can be considerable even in first mode when 
shaft and disk mass are comparable. 
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To study the effect of flexibility of disks, two 
shafts, one with single disk and other with two disks, shown in 
Table 6 are studied . Disk is treated as a standard shaft 
element. For the first shaft there are eight elements when disk 
is taken as rigid and nine elements when disk is flexible. 

For second shaft with two disks, there are eight elements when 
disks are taken as rigid and ten elements when disks are taken 
as flexible. 

The results are presented in Table 6. The critical 
speeds obtained assuming the disks to be flexible are higher 
than those obtained considering the disks to be rigid, as 
expected. In the first mode the difference in results is 5.2 */. 
for shaft with one disk and 8.1 '/.for shaft with two disks. 

3.5 DUAL SHAFT : 

A steel duel shaft with following parameters has been 
chosen for this study. 

0= =1.5, Re^= 0.05, C d (= -^) = 1.0 

2 i 

Critical speeds for this dual shaft with no shear def orimtion 
are obtained from Eqns . (2.55), and with shear deformation 
from Lqns . (2.53). For without shear case, results for forward 
and backward precession are obtained and are given in Table 7. 
With shear case results, for forward precession are obtained 
and are also given in the same table. 


Mode 


One disk 


Two disks 



Rigid 

Flexible 

Rigid 

Flexible 

1 

473.1 

488.6 

420.9 

457.8 

2 

1806.9 

1910.4 

1800.5 

1911 .1 

3 

4392.1 

4498.3 

3883.8 

4231 .7 

4 

8326.1 

8508.6 

8304.6 

8470.8 

5 

12357.2 

12824.8 

11506.6 

12392.5 




For forward precession, close forir solutions for 
ringed -hinged shafts are av..il~bie [4] end are given in 
Table 7, Present FSM results match ’..'ell with these. It is seen 
that dual stiffness has doubled the number of critical speeds 
when compared to the corresponding isotropic shaft. Further, 
as cn he seen from the eigen v lue problem Eqn. (3.9), for 
forv.'ard mode of procession, the- vibration in the two principal 
planes (Id and IIS) become uncoupled. Since the elastic 
properties are different in the two principal planes, there 
will be a critical speed critical- -speed corresponding to the 
el; stic property in each plane, loading to doubling of critical 
speeds in comparison to isotropic shaft. 


The critical speeds in backward mode of precession are 
also presented in Table 7. Also presented are, the critical 
speeds from [4], for hinged -hinged support conditions. By 
observing the eigen value problem, -oon. (3,12), it is seen 
that the vibrations in two planes are coupled by the [C-3 
matrix. Multiple critical speeds are seen to ex_st in 
backward mode, in contrast to forward mode, 

3.6 ISOTROMC SHAFT IITH DIFF^ i uIlIL 1 BUPrCKT CLafolT I OJ o 
TITTWO PLANTS : 

In this section, the results are obtained for cn 
isotropic shaft with different sup: ort conditions in two planes, 
delation is obtained from equations of dual shaft, nqn. (2.53), 
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lOule 7 : Critical speeds of a du~l shaft 
( 6-elements) 


without shear deformation (forward pression) 
! I 

^ ( c = = 1 .5 , Re^= 0.05, C d = 1.0) 


Mode 

Hinged -hinged 

Fixed -fixed 

Cantilever 

Present 

Exact 4 

1 

6,. 16 

S.16 

18.57 

2.89 

o 

it. 

y.99 

9.99 

22.73 

3.54 

3 

33.96 

33.95 

53.56 

18.77 

A 

41 .60 

41 .59 

65.60 

22.99 

5 

82 . 34 

82.30 

114.03 

56.11 

6 

100.64 

100.71 

139.66 

68.72 

7 

166.41 

165.74 

217.36 

123,18 

8 

203.81 

202.99 

266 . 20 

150.87 


With shear deformation 

(forward precession) 


(Me 

I 1 

— 2 .9 > T“ = 1.5 
x 2 

, Re/|= 0.05, 

= 1.0) 

Mode 

Hinged -hinged 

Present Exact \^~ 

Fixed -fixed 

Cantilever 

1 

7.96 

7.88 

16.57 

2.86 

2 

9 .6A 

9.54 

19.34 

3.48 

3 

30.61 

29 .44 

41 .87 

17.31 

A 

35.86 

34.55 

47.04 

20.45 

3 

64.07 

59.82 

74.61 

40.70 

6 

72.14 

67.74 

81 .35 

52.00 

7 

103.22 

94 ,33 

111 ,50 

82.18 

8 

112.30 

103.79 

118.88 

90.20 
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us ~ ]J a ~ 1. This sir: ft has Li 

v..uJ fired supports in IIw plane, 
conditions the forward critical s 
liven in Table £ . Also given are 
similar fixing conditions in two 


-'.ged sur.'Orts in IT plane 
U sing these b ound : ry 
peods are obtained and are 
the critical speeds, for 
planes i.e. similar fixing 


conditions of hin- cd 


and fixed 


sup sorts 


It is noted that critical speeds of different 
support conditions, HF-HF, can be obtained by combining the 
critical speeds of shafts of isotropic support conditions of 
IIH-IIII and FF-FF. Thus, there are two critical speeds in 
comparison to a shaft having simile r fixing conditions in the 
two planes of deflections. Critical speeds due to vibrations in 
each plane correspond to the support conditions in the 
respective planes. He ice it may be said that an isotropic 
shaft having dissimilar support conditions, behaves in a 
manner similar to a dual shaft having similar support 
conditions, with regards to critical speeds. 




»G e 

HH-HH 

HF-HF 

FF-FF 

1 

9 .64 

9.64 

19.34 

2 

35. 88 

1 9 . 9k 

47.64 

3 

72.14 

35.88 

81 .35 

4 

112.30 

47.04 

118.88 

5 

153.21 

72.14 

157.80 

6 

195.89 

81.35 

197.16 

7 

234.20 

112.30 

236.59 

C 

275.21 

118.86 

277.08 







65 


-' ; * 7 I^-^O-IC 31: .FT wTTII L.'.j TING : 

lo study the effect of internet viscous and 

hycleroLic damping, the Sons. (2.67) are solved. A steel 

chart with lie = 0.02 (i.e, i = 12.5) studied by [9} , with 

hinged -hinged support conditions is chosen. The bearings 

u ri i s: i med to have uncoupled stiffness and dumping 

charoctc-ristics i.e. If = r = 0 and C = C =0. 

hey yx xy yx 

dlir' complex frequencies are obtained in the form of 


^ .... (3,18) 

The nondimens ional frequencies ere plotted in 
. (3.1) to (3.5), for a speed range of 0 to 17 (i.e, 0 
to 1 Sod Orpin) , and for various combin'- tionc of, bearing 
flexibility and boarin 0 dumping characteristics. 

The curves labelled /, F”and * £ " represent forward 
rani backward precessional modes respectively. In all the 
casus studied, the results are compared with the results of 
[ 9 ] . The observations made in each case are listed 

below , 


Ons c 1 : Viscous dumping, isotropic undamped bearings. 

Bearings are undamped and have isotropic support 
conditions of 

K = K * 33.14 (17.5 MU/m) 
xx yy 

Considering the internal damping to be viscous, with 
= 0.01621 (0. 0002s), natural frequencies obtained are 
Jiven in Fig. (3.1) as function of speed q . The 
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Fig* 3.1 Natural frequency Aj vs. spin speed A 
for case 1. 





SI 


Fig-3-2 Natural frequency Aj vs. spin speed xl 
for case 2 
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Results from Nelson £9D 
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Fig. 3.3 Natural frequency A x vs. spih speed .n 
for case 3 
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Fig. 3-4 Natural frequency Aj vs spin speed A 
for case 4 
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iirr-L critical speed occurs at 6.44 (4980 rrra) and second 
crit.ic -1 speed occurs at 13.55 (10500 rpm) in forward mode 
ol precision, -^hese values are in good agreement with those 
given by [ 9 ] , which are 5.40 (4950 rpm) and 13.55, 

Iho first mode was found to become unstable at the first 
j. onward critical speed and the second mode at second 
lorv/rtrd critical speed. These modes stay unstable upto the 
ntu lied speed range (15.5). 


Co.se 2 : Viscous damping, isotropic damped bearings. 

The system studied in case 1 is analysed incorpor .ting 
isotropic bearing dampings with C , = C = 0.269 

XX vv 


(i.e. 1.75 x ~^~) where C 


yy 

represent the bearing 


m ' " xx’ yy 

dumping coefficients in non-dimensional form in X and Y 

directions respectively. The results are plotted in Fig. (3. 2). 

It is seen that, the first moce in now stable upto speed 

of 11 .6(9000 rpm) and the second mode is stable for the 

entire spin speed range studied (15.5). The results match 

very well with [ 9 }. ^ 

Case 3 : Viscous damping, anisotropic undamped bearings. 

Shaft has viscous damping = 0.0162. Bearings are 
undamped and support conditions ere anisotropic i.e. 


K 


xx 


33.140 , K. 


yy 


= 66.318, 
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The re:. alts are shown in Fig . (3.3). It is seen th: t the 
system becomes stable for the entire s_ in speed range 
studied. Forward and backward precession-:! speeds are wide 
. 1 ] art now. Results match very well with [9 ] 

C _ e 4 : Hysteretic damping, isomropic undamped bearing . 

The effect of hysteretic dnnroin_ is studied in this 
c.r.c. he' rin^j are undamped and have isotropic support 
conditions 

K ~ K = 33.14 
xx yy 

The value of loss factor r,^ used is 0.0002. The results are 
plotted in Fig. (3.4) . It is seen that both forward 
processional modes become unstable and both backward mooes 
arc stable for the entire spin speed range studied. As 
be Core forward and backward precessional speeds are quite 
close to each other in the modes studied. Results match 

very well with [9 ] • 

Case 5 : Hysteretic damping, isotropic damped bearings . 

Isotropic damped bearings with following properties 

arc used, 

K = K = 33.14 

K xx yy 

c = c = 0.269 

xx yy 


Value of loss factor ti used is 0.C002. 

The n. ul to ore shown in ^ig . (5.5). It is seen xhat both the 
modes ore stable for the entire speed range studied . The 

results match well witn [9 ] . Forward 'and backward preeesslonal speeds 
arc puite close to each other. 


Internal viscous damping is seen to have destabilising 
effect beyond the first critical speed in forward mode. The 
r, Lability of the system is improved by providing damped 
bearings. Stability is improved further by the provision of 
bu rings having anisotropic stiffness properties . Internal 
hystcretic damping destabilises all the forward moaes over the 
(■ntino speed ran^e . The provision of damped bearings is again 
soon to improve the stability of the system. 

5.7.1 If feet of bearing stiffness on stability, : 

For the shaft studied in section ;3.7),case 1, with the 
inotropic bearings, the instability was found to occure at the 
first critical speed. This cose was further extended to analyse 
the effect of bearing stiffness on the instability threshold, 
instability threshold speeds are obtained for a range of values 
of isotropic bearing stiffness. Tne results ^re 0 iven m 
Fig. (5.6). It is seen that, the instability threshold increases 
with increase in bearing stiffness significantly in the begihing. 


3.7 


7 Effect of bearing, damping on suability . 


♦ f • 

The shaft studied in Section (3.7.1), is analysed to 
study the effect of variation in bearing, damping properties at 
r few chosen bearing stiffness values, on the instability 



Fig, 3,6 Variation in instability threshold 
with bearing stiffness 
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threshold spend, for the first node. The* results obv-insd 
:ir - plotted in (3.7). It is seen that with an increase 

m besrin 0 damping coefficient the instability threshold 
increases linearly. The instability threshold is seen to be 
decreasing with an increase in stiffness of bearing s f when 
bcorin__ dumping is present. Thus a lower v: lue of bo ring 
sti.'fncss, 30.0, and higher value of bearing damping, 

0.35 (2.28 x 10^ ~p) , are seen to give the hig v est stable 
speed range in the present case. But as can be seen from the 
big . (3.7,) for a lower value of bearing damping coefficient, 
a beurin 0 with higher stiffness gives better stability. 




CONCLUSIONS 


The following conclusions emerge out of the 


present 


FEM is easily applicable to more familiar differential 
equations of rotor dynamics. It is simple and yields 
reliable results. Satisfaction of boundary conditions 
of intermediate supports (multi -span rotors) is 
s traight forward . 

When shear deformation effect is considered, results 
from two simultaneous second order partial differential 
equations using cubic polynomial approximation are 
very accurate (there is no ill conditioning) . 
Geometrical and natural boundary conditions are 
satisfied exactly. Results obtained are accurate with, 
fewer elements . 

Neglecting shaft mass can lead to large errors. The 
mass of the shaft is taken care of in a routine way 

in FEM. 

For undamped isotropic bearings, the instability 
threshold increases with increase in bearing stufness 
significantly in the begining . 
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For lower value of bearing damning higher 
stiffness increases the stability, and this 
effect is reversed for higher values of bearing 
aemping . 
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-■-PI UI2DIX - 1 


DERIVATION OF KINETIC EIU13Y 


The shaft cross section has a roc.tion 9 X about 
X-axis, 9 about Y-axis, ©„ about S-exis, translation 

y ^ 

x .-’long X-axis and y along Y-axis, as mentioned in section 
(2.1). The- kinetic energy due to rotation is obtained 
as follows 

Consider a frame of reference abc attached to shaft 
cross section such that, c axis is normal xo the cross section; 
a,b axes arc in the x lsne of cross seexaen and mutually 
perpaidicular to each other and to c axis. Then the rotations 
of cross section relaxive to XY3 frame are also the rotations 
of frame abc relative to XYS . Initially, the axes of abc 
frame are considered to be coincident with -ATS frame. The 
final position of abc frume relative to XYS is determined by 
the following sequence of rotations. 

1) A rotation by © x about X-axis gives ^ 2 h 2 c 2 

2) A rotation by 9 about b 2 ~axis gives 

3) A rotation by 9„ about -axis gives abc 


defer in 


to Fig. ( A -la) the position of a 2 b 2 c 2 frame 


after the rotation © x relative to XYS irame is 
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Fig. A-1 (a) Rotation by 6 X about X axis 

(b) Rotation by 8 y about b 2 axis 

(c) Rotation by 6 Z about c, axis 

(d) Final position of abc frame 
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► 
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0 
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- 
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(A-1.1) 

c 2 



0 

— 

- slne x 

C 0S& 

Jv 


o 

Uw» — ■ 




■Similarly the position of frame a^b^c^ relative to 
c^b o c 2 > cfter the rotation © is, Fig, (A-lb), 


Ol 

^1 


cos 0^ 0 -sin©y 


r — 

a 2 


) - 

oio 

< 

b 2 ^ 

~1 


sin© 0 cos© 

u y A, 


c 2 


Flir ?. ly, the position of fra. ns abc relative to frame a^b^c^ , 


after the rotation © z is, Fig, ('--1c) 


a 


cosG„ sinQ„ 0 

z 


cl^j 

- 1 b 

► = 

-sin©„ cosG 0 

Z 


bi > 

u C - 


o 

o 

_A 


_ c 1„ 



L J 



(A-1.3) 


The- "Inal position of abc after the above rotations is shown 
in Fi<: , (A —lb) . The frame abc has the an pals r velocities 

fe r along -a 2 , see Fig. (1-4 alon l b 1 and K along c * 

The components of \ along the directions a,b and c can be 
obtained by the following transformation using Sqns . (A-1.2) 

and (A-1 .3) 
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n 


cos 9^ 

sin9_ 

« 


-sin9„ 

cosG„ 

X,D 


z 

z 



0 

0 


0 

0 

1 


cos9 v 

J 

c 

- s in9 y 


. — 

- e x 

0 

1 

0 

< 
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f in6 y 

0 

cosGy 


0 

_ 


-9 V cos9 pos© 
x y z 


< @ x cos9y5in& 2 


-8 x sin9 y 


(A — 1 .A) 


8iiJi.ll rly the components 
Uiil c obtained using 


of G y along the directions 
Ipn. (1-1 .3) as 


a,b 


# 

9 


cosG sinG„ 0 


0 

y» a 


z ^ 



l*y.b 

> - 

-sin9^ cosG 0 

Ca 

< 

9 

y 

c' 


0 o 1 


0 

y,c 




_ _ 


© sinG^ 

y ^ 

1 e y cos6z ’ 


(1 -1 .5) 


0 



84 


ihe .mauler velocity vector is alon^ c direction, thus 

.iO u” ns.Lorirr. ti on is needed. Hence the '.njjler velocities 

clone 8 ,b and c directions, due to the rotations G , G 

x y 

and G„ are obtained using Euns . (1-1.4) end (A-1 .5) as 


Q 0 “ ©v a + ° v -- + ' 

D b = ®x,b + ®y,b +9 z,b 
°c = K c + K c + ®z c 


e x cosG y cos& z +© sin& 2 


G cos© sinG„ + & c os© 
x y z y z 


- 9 X sine y + G z 


(A-1 .6) 


..a the directions a,b and c are also the principal directions, 
tin Line tic energy due to trnnsl"tion rotation becomes 


T e = \ 7 [ m(x 2 +y 2 ) +I D (.of+ B 2 ) + I n o 2 ] ds (a- 1. 7) 
o 

Substituting Eqn. (A-1 .6), the expression for kinetic 
energy becomes. 

T= i J [m(x 2 +y 2 )+ I D (e 2 cos 2 G y + G 2 ) 
o ' 

+ I (G 2 sin 2 G +G 2 - 2G V G sin©,,) ] ds (A-1 .8) 
V x y ^ ^ j 


suming Q , 9 to be smell fro. no^in^ a. -t G^, ^ 


’x 7 y 


si In sneed of shaft), the above expression simplifies to 


% \ I [ m(X 2 +y 2 )+ l D (e 2 + e 2 ) 


2 , A2 


+ h ( Q o - 2 -°o ®x V l d5 (A' 1 -9) 
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a prai dp: - u 


DERIVATION OF DIFFERENTIAL ElUATIC. 


SY 1C., 'TOMS I STAGS 


Free body diagram of a differential element of shaft 
is shown in the Fig. (2.1). Newton’s law gives 


2F X = ma x $ (Q x + dQ^- Q x • mds x 

IHy = Icy (H y+ dM y )-H y + (V dQ x> Tf + Q x T- 

m It O d ds + J p ds D o 6 :A 


H? y = ma y ; (Q y «kj y )-Q y = mds y 

a ‘ x = C ; ('V dn x>- V V d V F ^ dE 


(A-2.2 ) 


d_ 

dt 



© +1 
x p 


dS C 0 



In these equations M represents "bending moment 
and Ci represents shear forces. Suffixes have the usual 
engineering sense in the right hand coordinate system, other 
terms have been explained in Chapter 2. angular momentums 
are shown in Fig.(A-2.) 

Neglecting higher order terms and rearranging, 
hens, (a -23) and (A-^2)become 
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'-C = mx 


K + U x 


I’-N ©_. + I Q © 

o y p o x 


(a-2.3) 


U y - my 


^ = “ Ip, + I 0^ ©. 


* y 


D x p o y 


U- 2.4) 


From mechanics of solids, it is known 


Q x - K ' GA c|f - v 

Q y ■ K ' GA < H - e x ) 


(A-2.5) 


K x = - EI 4 i f v = + EI e ) 


(A-2.6) 


Usinx these relations, Eons. (A*23) and (n-2,4)become 


[ k' ga (ft - ej]’ - mx = 0 


(El 0y)' + K'GA(|f - By) - I D S v -I p B 0 6 X = 0 


(A-2.7) 


[I'.'GA (|2 - 6 )] - my = 0 


(.El 6. 


dS X' 

k'ga(| 5 - ©J- i D e _+!_ o n e„ = 0 


(A -2. 8) 


X 


■3s x' jj x p “o y 


Equations (A -2.7) to (~ -2,8) are exactly the equations 
(2.9) and (2,10) obtained by Hamilton’s principle. 
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AFFE.'DIX- HI 


ST.-. TIC COKE 


-TION 


-hen the mass of a system is assumed to be concentrated 
-*t one or few chosen points, the mass matrix will be diagonal. 
Furtner, it contains zero diagorial elements for those 
degrees of freedom, for which there is no associated mass. 
V.'hcreas the element stiffness matrix formulated by FEM approach 
is fully populated and relates all the forces to corresponding 
degrees of freedom. But for those degrees of freedom, for 
which, the mass matrix contains zero element, the associated 
inertia forces are zero, hence, it is necessary to exclude 
those degrees of freedom from stiffness matrix, for dynamic 
analysis. The process is called as static condensation 
( Cough and Penzein [ 21 ]) and the procedure is as follows. 


Consider a typical matrix equation, 



# "Z. 

- - *” J 


D 


where represents the column of degrees of freedom to be 

condensed because the corresponding inertia forces are zero. 

represents the inertia forces associated with the 
degrees of freedom 
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From ^n. (n-3.1) 

[ K xJ W + [Kxy] “ {FjJ 

(A- 3.2) 

[Ky X ] ^ + [Kyy] \v] = < 0 > (A- 3 . 3 ) 

From Ain. (A-3.3) 

W = - tv 3 ' 1 C V 3 * x i (A - 3 - 4) 

Substituting this in 3qn. (A -3. 2) 

< [Kxxl - [A 3 Pyy 3 ’ 1 PyS > W ■ * F x5 

(At 3 . 5 ) 

Thus the degrees of freedom are condensed and expressed 

in terms of ' {x} and the condensed stiffness matrix [K Q ] 
is obtained as 

[K C ] = ([ K«] - [K^Pyy ]' 1 [r. yx ] ) (- 3 . 6 ) 
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